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STUDIES  OF  VON  KARMAN'S  SIMILARITY  THEORY 
AND  ITS  EXTENSION  TO  COMPRESSIBLE  FLOWS 

A  SIMILARITY  THEORY  FOR  TURBULENT  BOUNDARY  LAYER 
OVER  A  FLAT  PLATE  IN  COMPRESSIBLE  FLOW 
By  C.  C.  Lin  and  S.  F  .  Shen 

SUMMARY 


As  an  application  of  the  concepts  in  NACA  TN  2541,  the  prohlem  of 
turbulent  boundary  layer  over  a  flat  plate  in  compressible  flow  is 
treated.  ^  The  dissipation  term  in  the  energy  equation,  often  neglected, 
is  first  carefully  studied  and  found  to  be  of  importance.  In  the  dif¬ 
ferential  equations  governing  the  fluctuations,  the  lover- frequency  com¬ 
ponents  of  all  quantities  are  regarded  to  participate  in  an  equilibrium 
in  accordance  with  the  similarity  concept.  After  proper  linearization,  a 
set  of  differential  equations  containing  only  the  lower- order- frequency 
fluctuations  is  obtained.  In  parallel  with  Von  Karman’s  theory  in  incom¬ 
pressible  flow,  the  similarity  scales  for  all  the  flow  variables  are 
derived.  Two  possible  length  scales  are  found,  and  the  significance  of 
this  possibility  is  discussed. 


INTRODUCTION 


Prior  to  the  development  of  high-speed  aircraft,  the  aerodynamic 
phenomena  were  satisfactorily  explained  by  theories  based  on  an  incom¬ 
pressible  fluid.  Such  a  situation  certainly  no  longer  remains  true. 

An  increasing  abundance  of  research  is  now  available,  attempting  to 
extend  the  theories  into  the  compressible  range  for  almost  every  aspect 
of  practical  importance.  Yet  not  until  recently,  has  the  problem  of 
turbulent  boundary  layer  in  compressible  flow  attracted  the  attention 
of  investigators.  The  reason  is,  of  course,  that  any  turbulence  problem 
is  a  difficult  one  even  without  compressibility.  ,  Perhaps  demanded  by  its 
practical  importance  in  the  flows  over  an  airfoil  or  in  a  wind  txmnel, 
or  just  to  keep  pace  with  the  advancement  of  knowledge  in  other  aero¬ 
dynamic  problems,  papers  dealing  with  the  behavior  of  turbulent  boundary 
layer  in  compressible  flow  finally  begin  to  appear. 
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Thus  one  has,  for  instance,  the  works  of  Ferrari  (reference  l) , 
Wilson  (reference  2),  Van  Driest  (reference  3)^  Ladenburg  and  Bershader 
(reference  4),  and  Eckert  (reference  5)*  Except  for  reference ^4,  which 
is  a  pure  experimental  measurement  using  interferometric  technique,  all 
of  them  present  a  theoretical  analysis.  On  reviewing  these,  it  appears 
that,  in  these  analyses,  the  turbulent-dissipation  tenns  were  neglected 
by  the  argument  that  any  molecular  phenomenon  is  probably  of  no  impor¬ 
tance  in  comparison  with  the  turbulent  transfer.  This  argument,  being 
without  experimental  backing,  should  be  carefully  investigated. 

The  first  part  of  this  report  is  concerned  with  such  an  investiga¬ 
tion.  The  energy  equation  is  studied  in  some  detail,  and  the  orders  of 
magnitude  of  the  various  terms  are  estimated,  with  the  help  of  available 
experiments.  The  turbulent  dissipation  is  found  to  be  not  negligible  at 
least  for  boundary- layer  flow  over  an  insulated  wall.  Von  Karman*s 
formula  for  the  mixture  length  has  been  used  in  this  analysis;  but  the 
conclusions  reached  on  the  order  of  magnitude  of  the  various  terms  are 
presumably  correct,  even  though  the  distribution  of  scale  of  turbulence 
may  be  only  roiighly  given  by  that  formula. 

The  second  part  of  this  report  is  an  attempt  to  extend  Von^Karman’s 
similarity  theory  to  the  turbulent  flow  of  a  gas  in  the  boundary  layer. 
This  extension  is  based  on  the  ideas  developed  in  reference  6.  It  turns 
out  ( see  reference  7)  that  the  theory  has  a  certain  universal  character 
with  regard  to  the  Mach  niomber  and  the  heat-transfer  conditions  at  the 
wall,  and  the  unavoidable  correlation  constants  can  be  determined  once 
for  all.  Owing  to  the  lack  of  suitable  experimental  data,  the  numerical 
part  of  the  theory  cannot  be  pursued  and  is  left  to  the  future. 

This  work  was  conducted  at  the  Massachusetts  Institute  of  Technology 
under  the  sponsorship  and  with  the  financial  assistance  of  the  National 
Advisory  Committee  for  Aeronautics. 


ENERGY  EQUATION  AND  DISSIPATION  TERM 


In  problems  of  compressible  flow,  a  relation  between  the  velocity 
field  and  the  temperature  field  is  furnished  by  the  equation  of  energy. 

In  applying  the  energy  equation  to  turbulent  phenomena,  it  has  been  the 
customary  practice  to  neglect  both  the  heat  conductivity  and  the  vis¬ 
cosity  terms.  The  reason  implied,  or  sometimes  stated,  is  that  the 
molecular  phenomena  would  make  but  very  small  contributions  in  comparison 
with  the  turbulent  transfer  terms,  just  like  the  situation  in  the 
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equations  of  momentum.  Undoubtedly  such  vague  arguments  have  definite 
limitations.  The  earlier  applications  were  for  the  problem  of  tempera¬ 
ture  distribution  in  turbulent  flow  at  low  mean  velocities.  There  is  a 
predominant  temperature  field.  The  conduction  terms  are  of  the  same 
nature  as  the  viscous  shear  in  the  equations  of  momentum.  The  heat 
generated  through  viscous  dissipation,  being  caused  by  the  velocity 
fluctuations,  is  probably  small  because  of  the  large  mechanical  equiv¬ 
alent  of  heat.  Hence  it  seems  reasonable  that  for  low  Mach  numbers  the 
neglect  of  molecular  terms  would  be  justified.  Recently,  in  treating 
turbulent  boundary  layer  in  supersonic  flow,  various  authors  have  retained 
such  a  procedure  without  questioning  its  validity  (references  1  and  3)  • 

The  following  is  an  attempt  to  disclose  the  deciding  factors  in  assessing 
the  relative  weight  of  the  molecular  terms  and  the  turbulent  terms.  Most 
important  of  all,  the  dissipation  term  is  found  to  be  not  negligible  in 
the  case  of  boundary-layer  flow  over  an  insulated  wall,  irrespective  of 
Mach  number. 

One  may  first  write  down  the  complete  energy  equation  (reference  8, 

p.  6o6) 

DT  (bvL  .  8v  .  5w\  A. 

5t  *  ^  ®  (1) 


where  c^  is  the  specific  heat  at  constant  volume,  k  is  the  heat 
conductivity,  and  0  is  the  dissipation,  defined  by 


$  =  (1 


(2) 


The  first  term  on  the  left-hand  side  represents  the  convective  action. 
Equation  (l)  is  the  equivalent  of  the  statement  that  following  a  particle 
the  temperature  change  is  due  to  three  kinds  of  heat  sources:  The  com¬ 
pressibility,  the  heat  conduction,  and  the  dissipation.  Without  essential 
difference,  let  the  incompressible  case  be  considered.  Then  equation  (l) 
becomes. 


1.^  =  1^ 
7V  Dt  CT 


(3) 
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obtained  by  omitting  the  compressibility  terms  and  dividing  through  with 
Cp^L.  (See  appendix  for  definitions  of  symbols.)  Here  7  =  Cp/cy^  the 

ratio  of  specific  heats,  a  =  Cp^i/k,  the  Prandtl  number,  and  e  = 

In  turbulent  flow,  each  quantity  may  be  split  into  mean  and  fluctuating 
parts;  that  is, 


1  i 

/dT 

1  =  +  iv^T' 

+  -i-(6  +  e'] 

7v' 

^Dt  ^  Dt  j 

f  a  a 

(3a) 


It  is  now  intended  to  compare  the  terms  on  the  right-hand  side  with  a 

dT 

typical  term  on  the  left-hand  side,  such  as  v*  — .  The  mean  flow  is 

assumed  to  have  a  characteristic  velocity  U  in  the  x- direction,  as  in 
boundary-layer  or  wake  flows . 


(1)  The  molecular  conduction  term 

Let 


due  to  mean  temperature: 
o 


Then 


i-l  =  V 


°  dy7’'^ 


V  1  /d^/dT\ 
o  V '  \dy  dy  / 


(in  the  following  discussions  of  the  orders  of  magnitudes  suitable  root- 
mean-square  values  will  be  implied  without  explicit  notation.)  At  least 
for  a  study  of  the  order  of  magnitude,  one  may  use  Reynolds’  analogy 
(reference  8,  p.  649  ff>)  to  relate  mean  temperature  with  mean  velocity; 
that  is. 
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where  Z  may  be  interpreted  as  either  the  mixing  lei^gth  or  Von  Karman*s 
similarity  scale.  Thus^ 

(4a) 


or 


(2l  JLh  1.) 

[a  v'  Z  rJ 


m 


where  L  is  the  characteristic  length  of  the  mean  flow  and  Rj^ 
mean-flow  Reynolds  number,  defined  by 

Rl  =  im/v 

1 

Two  illustrative  cases  may  be  given.  For  jet  or  wake,  flow,  L  may  be 
taken  as  the  "breadth”  of  the  wake,  where  U  =  ~  The  well-known 

theories  (see,  e.g.,  reference  9)  assume 

L/Z  =  Constant 

the  constant  being  of  the  order  of  0.2  for  two-dimensional  jets.  The 
ratio  r2^  is  indeed  small  throughout  the  section.  In  the  case  of 
boundary-layer  flow,  L  may  be  taken  as  the  thickness  of  the  layer  and 

Z  0.4y 

Then  equation  (4b)  indicates  that  r]_  becomes  important  only  when  very 
close  to  the  wall,  where  laminar  sublayer  comes  in  anyway. 

A  different  form  of  formula  (4a)  will  also  show  that  r^^  is  gen¬ 
erally  small.  Introducing  the  turbulence  Reynolds  number  R^  =  u*A,/v, 
where  X  is  the  microscale,  equation  (4a)  becomes 
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(2)  The  molecular  conduction  term  iv^T'  due  to  fluctuating 

a 

temperature:  The  term  involves  second  derivatives  of  a  fluctu¬ 

ating  quantity,  therefore  greatly  influenced  by  the  high-frequency 
components  of  the  fluctuation.  It  will  be  compared  with  a 
1  ' 

term  —  v'  - .  Before  the  comparison,  however,  one  may  observe  that 

TV 

the  neglect  of  viscous  shear  in  the  equation  of  momentum  implies  the 
following  relation: 


^2 


«1 


(5) 


in  spite  of  the  high-frequency  components  in  S^u’/^y^-  form  the 

ratio 


r  =  Iv^T'/—  v' 

^  a  j  7V  oy 

^  V 

0  v'W^/ay  j 


(6a) 


Rewriting, 


^3 


Su*/Sy  \ 
ST*/&y  S^u’/Sy^/ 


(6b) 


Reynolds^  analogy  assximed  that  the  temperature  and  velocity  fluctuations 
are  proportional.  Although  such  an  approximate  picture  cannot  be  expected 
to  hold  for  the  derivatives  largely  controlled  by  the  high-frequency  com¬ 
ponents,  it  is  likely  that,  at  least, 

A  ^/i  ^ 

AT  ay2  /  U  ay2 

4|r/l|u:,0(l) 

AT  9y  /  U  5y 
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where  AT  is  the  mean  temperature  difference  across  the  boundary  layer 
and  U  is  the  mean  free-stream  velocity.  Then  equation  (6b)  reduces  to 

rj  .  o(^),  «  1 


(3)  The  dissipation  terms  "?  and  The  order  of  magnitude  of 

the  dissipation  term  is  to  be  compared  with*  v*  dT/dy  on  the  left-hand 
side.  To  begin  with^  one  discards  the  contribution  of  the  mean  flow 
velocity  gradients  to  the  dissipation,  as  a  consequence  of  the  discussion 
in  reference  6.  Then, 


The  relative  magnitude  of  and  e  may  be  seen  by  studying  the 

following  equation: 


for  both  isotropic  turbulence  and  in  the  wake  of  a  cylinder.  Since  the 
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high-frequency  components  are  responsible  for  €,  one  expects  the  ratios 
in  the  two  cases  to  be  close  to  each  other.  It  turns  out  that  the  wake 
case  has  a  higher  ratio  but  both  were  within  10  percent  of  the  value 
It  may  therefore  be  concluded  that 

0(l) 


For  the  order  of  magnitude  of  the  mean  dissipation  €  Taylor's 
work  may  first  be  cited  (reference  12).  For  a  channel  flow  and 
omitting  compressibility,  the  work  done  on  an  elementary  volume  is 

.T  — .  where  t  is  the  total  shearing  stress,  viscous  plus  turbulent, 
dy*^ 

Based  on  experimental  measurements,  the  microscale  was  evaluated  and  the 
dissipation  estimated  by  the  formula 


e  15i-L 


(u')g 


which  was  originally  derived  for  isotropic  turbulence.  Plotted  across 
the  section,  his  results  indicate  that  these  two  terms  are  essentially 
of  the  same  order  of  magnitude  over  a  greater  part  of  the  section  except 
near  the  center  line  and  the  wall.  The  dissipation  further  does  not 
change  its  order  of  magnitude  in  these  regions.  The  transfer  of  turbulent 
energy,  being  the  difference  of  the  two,  is  therefore,  in  general,  not 
sufficiently  important  to  alter  their  orders  of  magnitude.  It  seems  now 
justifiable,  based  on  Taylor's  experience,  to  estimate  the  order  of 

magnitude  of  7  by  evaluating  ^  points  near  half  the  thickness 

in  the  case  of  a  boundary -layer  flow,  where  the  velocity  profile  is 
rather  similar  to  a  channel  flow.  In  a  jet  flow,  the  region  of  maximum 
viscous  shear  might  be  chosen  for  the  estimation  in  order  to  avoid  the 
nearly  isotropic  turbulence  in  regions  of  weak  shear.  The  correspondence 
there  is  undoubtedly  less  satisfactory  but  perhaps  still  close  enough 
for  the  present  purpose. 


It  also  may  be  easily  shown  that,  in  the  case  of  axial  symmetry, 
the  work  done  per  unit  volume  becomes  t  Taking  this  to  represent 

_  ,  St 

the  order  of  magnitude  of  e,  one  may  compare  with  the  term  pc^Up  ^ 
to  determine  the  relative  importance  of  e  in  the  energy  equation. 


^For  a  Gaussian  distribution  of 


3.0.  Experiments  by  Simmons  and  Salter  (reference  ll) 
follows  very  closely  a  Gaussian  law  and  indeed  gave  (u') 


Su'/Sx,  the  ratio  would  have  been 
indicated  that  u 
2 


7[{u’)2] 


3.0 
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To  estimate  t,  when  the  mean  motion  depends  essentially  on  y  and 
the  normal  velocity  v  o,  , 


-r  .  du  - 

^  ^  puv 


du 


-^  +  p  u'v' 
dy 


for  incompressible  flow.  Except  very  close  to  the  wall,  t  p  u'v'. 
This  formiaa  may  be  regarded  as  valid  even  in  compressible  flow  if  only 
the  order  of  magnitude  is  desired.  Rewriting, 


where  is  the  correlation  function.  To  compare  the  dissipation 

with  transfer  terms  in  the  equation  of  motion,  form  the  ratio 


"•5  = 


dy 


Then, 


r^  Rs 


du 


Ruv 

u* 

Um^ 

Urn 

T™ 


(9a) 


where  Ujjj  is  the  characteristic  mean  velocity,  Tjjj  is  the  characteristic 

mean  temperature,  and  ATj)j  is  the  characteristic  mean  temperature  differ- 

du  /dT  _  _ 

ence.  The  replacement  of  — /—  by  Um/ATm  is  again  based  on  Reynolds’ 


10 


MCA  TN  25*^2 


analogy,  identifying  the  mean  velocity  and  temperature  distributions. 
With  the  introduction  of  a  characteristic  Mach  number  Mjjj  equation  (9a) 
becomes 


(9b) 


Thus  it  is  seen  that,  with  a  similar  turbulence  pattern,  the  importance 


(AT. 


m 


of  the  dissipation  term  is  proportional  to  /-zr—*  The  similarity  of 


■m 


the  turbulence  pattern,  incidentally,  needs  only  to  cover  the  lower  part 
of  the  frequency  spectrum,  as  would  prevail  when  mean  streams  of  differ¬ 
ent  velocities  are  passed  through  the  same  grid. 


The  criterion  (9b)  will  now  be  applied  to  two  cases  of  experiments  on 
jet  flow  made  by  Corrsin  (references  13  and  l4) ,  In  reference  13,  the 
data  are  as  follows: 


Jet,  1  inch 

Uq  =  Nozzle  velocity  =  10  meters  per  second 
ATq  =  Nozzle  temperature  difference  ^  10^  C 

At  x/d  =  20,  r  =  5  centimeters  from  axis 
^m/ ~  ^ 

^m/ATo  =  0.22 

|u'/%|  0.l8 

|v’/Ujii|  0-1^ 
u'v'/um^  =  0.010 

^(u/Um)  0.132  per  centimeter 


0.120  per  centimeter 
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(The  slopes  are  almost  identical,  indicating  the  approximate  validity 
of  the  Reynolds'  analogy  in  this  case.)  Hence, 


But 


Thus 


rpr  0  0.07  — 


yim 


6k  X  10*^ 
^m/Tm  T  X  10-3 


0(0.001) 


In  reference  l4  hy  Corrsin  and  Uberoi  there  are  the  following  data 

Jet,  1  inch 

Uq  =  100  feet  per  second 

ATq  =  170°  c 

At  x/d  =  15,  ATjh  0.3 
aTq  =  50°  C 
Ujn  O.k 

Uq  =:  4o  feet  per  second 

At  t/t2  =  1.0  of  the  section,  where  r2  is  the  radius  at  which 
^/ATja  =  1/2, 

^(u/Hn)  <«  0.6  per  centimeter 

^^AT/aTih)  «<0.55  per  centimeter 

|v'/i^|  ^0.1 

u'v'/i^2  «« 0.006 


I 
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Hence 


rc  0(0.06  -1.  _ 

ATjn/Tm 


With  0.0016  and  ATju/Tm  0.l4,  one  has  finally,  again, 

0(0.001).  The  experience  in  these  two  cases  seems  to  suggest  that 
in  jet  flow  the  dissipation  is  of  no  importance  in  the  energy  equation. 

For  qualitative  results,  equation  (9fc)  may  he  put  into  another  form. 
Introducing  the  mixing  length, 

u^  =  llv'lfi 

dy 

Then 


rc-  0 

5 


Z 


d  /  u  \ 

\Um/ ^m/”l^m 


(10)  , 


This  formula  may  he  checked  hy  the  Jet-flow  results.  If  the  "hreadth"  h 
is  taken  to  characterize  the  ^t  spreading,  similarity  of  the  mean  veloc¬ 
ity  profile  will  lead  to  =  A/h  at  any  given  y/h,  say  y/h  =  1. 

Using  Prandtl's  assumption  of  2  in  jet  flow, 

2  =  Bh 

Consequently,  eq-uatlon  (lO)  becomes 


r^  A  X  B  X 


(11) 


From  previous  experiments,  A  X  B  0.06j  its  constancy  is  thus  verified. 

Let  equation  (lO)  be  now  applied  to  boundary -layer  flow.  Here 
may  take 


Mjn  =  Free -stream  Mach  number 
Tjjj  =  Free -stream  temperature 


one 
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Consider  the  case  of  an  insulated  wall, 

ATjji  =  Stagnation  temperature  -  Tm 


as  usually  accepted.  It  follows  that 

y_  ~  g 

AVTm  7  -  1 

a  constant,  and 


(12a) 


Now  in  incompressible  flow^  a  popular  approximate  form  of  the  velocity 
distribution  is 


Urn  U/ 


where  6  is  the  boundary -layer  thickness.  Ladenburg  and  Bershader 
measured  for  a  supersonic  stream  over  a  flat  plate  (reference  ^4-), 
concluding 


Putting  u/i^  =  (y/5)^,  and  taking  Z  =  Ky  as  for  incompressible  flow, 
one  gets  for  order -of -magnitude  purposes 


(12b) 
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Ab  previously  stated,  the  dissipation  may  he  estimated  at  the  station 
y/6  0.5*  With  K  OA,  7  =  1.4  for  air, 

r^  ^  o[2n(0.5)^ 

0(0.3)  for  n  =  1/7 

-  0(0.25)  for  n  =  1/9  (13) 


The' result  of  equation  (13)  shows  that  for  an  insulated  wall,  the  dis¬ 
sipation  is  of  the  same  order  of  magnitude  as  the  terms  kept  in  the  energy 
equation.  The  usual  practice  of  omitting  the  dissipation  term  is  there¬ 
fore  incorrect.  It  is  also  interesting  to  see  that,  in  boundary -layer 
flow  over  an  insulated  wall,  the  weight  of' the  dissipation  term  is 
sensibly  independent  of  the  Mach  number.  Such  a  situation  naturally  will 
be  greatly  modified  when  excessive  cooling  or  heating  of  the  wall  is 
Introduced  by  auxiliary  means. 


EXTENSION  OF  SIMILARITY  THEORY  TO  COMPRESSIBLE 
TURBULENT  BOUNDARY  LAYER 


After  the  exploration  in  reference  6  of  the  foundations  and  limita¬ 
tions  of  the  similarity  concept  in  the  case  of  incompressible  turbulent 
boundary-layer  flow,  it  is  natural  to  attempt  to  formulate  a  theory  of 
the  compressible  turbulent  boundary  layer  on  a  somewhat  similar  basis. 

One  recognizes  that  the  additional  density  and  temperature  fields,  the 
mean  distribution  as  well  as  the  fluctuations,  must  bring  in  many  more 
difficulties.  The  similarity  theory  was  shown  to  involve  many  approxi¬ 
mations  for  the  incompressible  problem.  The  extension  to  compressible 
flow  must  therefore  be  made  to  an  even  more  approximate  degree.  Never¬ 
theless,  it  may  again  be  stressed  that  the  similarity  concept  leads  to 
a  simple  model  of  the  intrinsic  turbulence  mechanism,  which  is  essentially 
correct  at  least  in  incompressible  flow.  By  using  this  model  a  unified 
theory  is  possible  for  the  mean  velocity  and  temperature  phenomena,  with¬ 
out  the  necessity  of  separate  assumptions  whose  consistency  with  each 
other  cannot  be  ascertained.  Moreover,  the  mean  distributions  are  not 
very  sensitive  to  the  assumed  mechanism  of  the  turbulence,  as  the  exper¬ 
ience  in  compressible  flow  has  indicated.  Any  discrepancy  of  the  simi¬ 
larity  model  from  the  true  one  may  be  expected  to  Influence  but  little 
the  usefulness  of  the  theory  in  predicting  mean  distributions. 
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A  review  of  prevailing  theories  of  compressible  turbulent  boundary 
layer  discloses  that  a  better  fundamental  concept  can  be  claimed  by 
none.  The  earliest  attempt  by  Von  Karman  (reference  15)  was  to  sub¬ 
stitute  the  wall  conditions  for  the  stream  conditions  in  the  formulas 
derived  for  incompressible  flow.  More  recent  ones  include  the  compar¬ 
atively  simple  theories  of  Wilson  (reference  2)  and  Van  Driest  (refer¬ 
ence  3) ^  both  starting  from  the  incompressible  expression 


length  Z  =  Ky,  while  Wilson 


formally  carried  out  by 
allowing  p  to  vary  as  a  consequence  to  the  temperature  distribution. 

By  assuming  the  effective  turbulent  Prandtl  number  equal  to  unity,  the 
mean  temperature  and  the  mean  velocity,  in  the  case  of  an  insulated  wall, 
satisfy  the  isoenergetic  law, 

— •  1  •— o 

CpT  +  -  u^  =  Constant 

in  an  analogous  manner  as  for  laminar  flow.  A  recovery  factor  (refer¬ 
ences  16  and  17)  was,  sometimes  employed  as  a  refinement  of  the  Isoener¬ 
getic  law.  For  empirical  data  on  the  effective  turbulent  Prandtl  number 
and  the  validity  of  the  isoenergetic  law  in  turbulent  flow,  discussions 
can  be  found  in  literature.  In  Jet  flow,  for  instance,  the  effective 
Prandtl  number  for  air  lies  between  O.7  and  0.8  (reference  18) .  But  if 
dissipation  terms  are  not  negligible,  the  isoenergetic  law  cannot  follow 
in  any  case.  Though  useful,  it  must  be  regarded  as  largely  empirical. 

The  work  of  Ferrari  (reference  l)  started  on  a  more  solid  basis  by 
including  the  energy  equation  to  furnish  a  relation  between  the  mean 
temperature  and  velocity.  The  assumption  of  a  special  relation,  such 
as  the  isoenergetic  law  for  the  insulated  wall,  is  thus  avoided.  He 
used  Prandtl ’s  mixture -length  concept  to  deal  with  momentum  and  enthalpy 
transfers.  When  the  two  mixture  lengths  are  taken  to  be  identical,  a 
temper at  lire -velocity  relation  is  obtained  and  appears  to  be  formally  the 
same  one  which  occurs  in  a  laminar  boundary  layer  for  a  Prandtl  number 
of  unity.  This  coincidence  perhaps  may  be  used  to  explain  the  fact  that 


T  =  pZ' 


du 

dy 


du 

dy 


Van  Driest  took  Z  to  be  Prandtl 's  mixing 


used  Von  K^man’  s  form 


,  «  du  /d^ 

dy/dy2 


The  extension  to  the  compressible  case  v;as 
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the  isoenergetic  law  has  been  found  to  lead  to  no  serious  error.  However^ 
Ferrari,  too,  neglected  the  viscosity  terms  in  both  the  equations  of 
momentum  and  that  of  energy.  His  procedure  cannot  be  easily  applied  if 
the  dissipation  terms,  which  are  shown  to  be  not  negligible,  are  kept. 

In  addition,  his  argument  for  identifying  the  two  mixture  lengths  (for 
momentum  and  for  enthalpy)  is  by  no  means  conclusive. 

It  is  thought  that  especially  in  the  compressible  case  an  intuitive 
guess  to  the  proper  form  of  the  mixture  length  is  rather  difficult  because 
of  the  complicated  interactions  of  the  velocity,  density,  and  temperature. 
A  systematic  approach  based  on  a  clearly  defined  mechanism  becomes  nec¬ 
essary  in  order  to  achieve  self-consistency. 

As  was  previously  shown  in  the  discussion  of  the  incoii5)resslble  case, 
a  theory  can  be  established  on  an  assumed  similarity  among  the  lower- 
frequency  components  of  the  turbulent  fluctuations.  In  extending  to  the 
compressible  case,  the  behaviors  of  the  density  and  temi)erature  fluctu¬ 
ations  require  certain  additional  simplifications  based  on  physical  con¬ 
cepts.  The  differential  equations  are  to  be  satisfied  in  a  more  approx¬ 
imate  manner.  Let  the  complete  set  of  equations  be  first  written  down: 


Dui 

Dt 


1 

—  T— ^  i 
P  hx± 


1,2,3 


m 


1  Dp  ,  I  I 
p  Dt  5x  Sy  5z 


(15) 


Sx  5y  Sz. 


= 


(16) 


p  =  pRT 


(IT) 


One  may  note  that  in  the  equations  of  momentum  the  viscosity  terms  are 
neglected  but  not  in  the  equation  of  energy  as  a  consequence  of  the 
discussion  in  the  section  "Energy  Equation  and  Dissipation  Term."  Prior 
to  a  mathematicfij.  derivation,  the  procedure  may  be  summarized  to  consist 
of  the  following  considerations: 

(l)  Partial  linearization  of  the  differential  equations  by  neglecting 
the  fluctuations  of  density  and  temperature  against  the  mean  quantities 
themselves  but  not  the  derivatives  of  the  fluctuations. 


I. 
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(2)  The  separation  of  the  density  and  temperature  fluctuations  into 
lower-  and  higher -frequency  components  in  a  manner  analogous  to  the 
velocity  fluctuations.  A  similarity  of  the  turbulence  is  then  assumed 
to  include  all  these  lower -frequency  conqponents  • 

(3)  The  dilatation  of  fluid  element,  that  is. 


e 


5u  ^  Sv  ^  5w 
hx  Sy  hz 


is  regarded  as  a  separate  entity,  not  deduclble  from  simil^ity  of  the 
lower-frequency  components. 


The  following  discussions  will  substantiate  these  points: 


The  partial  linearization  process  is  based  on  the  fact  that,  ordi¬ 
narily,  the  magnitude  of  turbulent  fluctuations  is  only  a  small  fraction 
of  the  mean  quantity.  For  a  first-order  approximation,  one  may  take 


r*  ”1 

wm 

,  «1 

► 

/(t')7t 

,  «1 

(18) 


In  takixig  derivatives,  however,  the  fluctuation  terms  must  be  kept  as 
was  previously  done  for  the  velocity  fluctuations  and  with  the  same 
arguments.  By  so  doing  it  is  seen  that  the  differential  equations  are 
satisfied  to  a  lesser  degree  in  comparison  with  the  incompressible  case. 
The  difference  between  velocity  and  other  quantities  lies  in  that  the 
mean  density  and  temperature  cannot  be  dismissed  like  the  mean  velocity 
by  the  introduction  of  relative  motion.  The  mean  density  distribution 
alone  essentially  provides  the  coupling  between  temperature  and  veloci¬ 
ties.  A  consequence  of  equations  (18)  is  that  the  equations  of  momentum 
will  look  exactly  like  those  in  incompressible  flow,  only  with  a  variable 
mean  density  distribution.  Thus,  as  very  rough  approximations,  the 
starting  point  of  the  theories  of  Wilson  (reference  2)  and  Van  Driest 
(reference  3)  is  not  without  justification. 
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For  a  closer  examination^  one  separates  the  larger  eddies  from  the 
small  ones. 


P'  =  Pz'  +  Ph' 
T'  =  T2'  +  Tji' 


(19) 


where  subscripts  Z  and  h  denote  the  lower-  and  the  higher -frequency 
components,  respectively.  The  fluctuations  are  then  seen  to  be  due  to 
two  sources,  convective  actions  of  the  larger  eddies  and  heating -from 
viscous  dissipation.  Loosely  speaking,  the  convective  actions,  being 
carried  by  eddies  of  the  larger  sizes,  would  show  up  only  in  the  lower- 
frequency  components.  The  viscous  heating  comes  from  the  very  small 
eddies  and  would  contribute  to  both  the  lower-  and  the  higher- frequency 
components.  The  demarcation  between  the  "larger"  and  "small"  eddies 
therefore  becomes  even  less  clean-cut  than  that  in  incompressible  flow, 
where  the  coupling  through  a  density  variation  does  not  exist.  It  is 
nevertheless  still  conceivable  that  for  eddies  small  enough,  their  equi¬ 
librium  is  controlled  by  no  other  than  Kolmogoroff * s  (reference  19) 
parameters  jl  and  e,  fl  being  now  the  local  mean  viscosity.  For  points 
across  the  boundary -layer  thickness,  there  is  merely  a  shift  of  the 
Kolmogoroff  range  in  the  wave -number  space,  the  ratio  of  the  wave  numbers 
separating  the  Kolmogoroff  range  from  the  rest  being  of  order  unity  since 
the  ratio  of  \I  is  only  of  order  unity.  Tlie  higher -frequency  components 
p^’  and  Tj^'  therefore  must  participate  in  a  different  equilibrium  from 

and  T^^.  The  attempt  now  is  to  establish  differential  equations 
governing  the  approximate  behavior  of  the  lower -frequency  component  only. 

It  is  found  that,  after  partial  linearization,  a  typical  term  to 
be  examined  in  the  process  of  setting  up  the  differential  equations  for 

the  lower-frequency  fluctuations  is  ^  T'u' .  Reasoning  in  the  same  way 

ox 

as  before,  the  lower -frequency  part  of  T’u'  consists  of  T^’u^’  plus 
some  contribution  from  T^'u^^'  .  Now  the  magnitude  of  Tj^’  is  probably 
of  the  same  order  as  T^r’.  For,  one  could  visualize  that  T^’  for  order 
of  magnitude  is  associated  with  the  mean  distribution  dT/dy,  while  Tj^' 

is  essentially  related  to  the  dissipation  c’.  The  quantity  e’  was 
shown  in  the  section  "Energy  Equation  and  Dissipation  Term"  to  be  of  the 
s^e  order  as  6  and  also  of  the  same  order  as  the  contribution  from 
dT/dy.  However,  previous  discussions  indicated  that  «  ^U2’^^j 

hence 


\'^h'  =  Tz’u^'  uj,' 


«  1  (20) 
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Therefore,  the  high-frequency  parts  again  may  he  omitted  in  terms  like 
T'u'  without  seriously  affecting  the  equilibrium  of  the  lower -frequency 
components. 


There  is  another  way  of  looking  upon  such  terms  as  ^  T'u’ 

ox 

occurring  in  the  energy  equation^  relying  more  upon  a  physical  inter¬ 
pretation.  One  could  regard  the  operator  D/Dt  as  representing  the 
change  following  the  fluid  element.  The  higher -frequency  velocity  fluc¬ 
tuations  tend  to  average  out  without  contributing  appreciably  to  the 
movement  of  the  element. ^  The  operator  D/Dt  should  therefore  be  asso¬ 
ciated  with  only  u^'^  ^  equation  (l^),  the  quantity 

after  D/Dt  is  the  velocity.  Being  now  lineax,  the  equation  can  be  ' 
split  into  one  for  the  higher  frequency  uj^'  and  another  for  u^*.  The 

previous  results  for  the  velocity  fluctuation  are  reproduced.  In  eq\ia- 
tion  (l6),  the  temperature  appears  also  in  a  linear  fashion,  the  lower- 
and  higher- frequency  parts  can  again  be  slB5>ly  separated. 


The  dilatation 


Su  5v  Sw  .  n 

^  ^  deceptively  simple  in  appearance^ 


w  - 

requires  a  different  handling  when  a  theoiy  of  the  lower- frequency  com¬ 
ponents  is  desired.  Firstly,  the  magnitude  of  each  term  in  e  is  mainly 
due  to  the  high-frequency  component  rather  than  the  low,  since 


0 

Y ' 


has  been  shown  in  reference  6  to  depend  on  the  very  small 


eddies  in  the  study  of  the  dissipation  e .  In  the  second  place,  the 


2 

One  may  get  an  idea  of  the  relative  importance  to  diffusion  of  the 
low-frequency  and  the  high-frequency  components  by  examining  Taylor's 
formula  (reference  20) 


i  y2  =  (v-  ^  ^  dt 

for  diffusion  by  continuous  movements  where  R|  is  the  correlation  of 

the  velocities  at  time  interval  I  apart.  Being  proportional  to  (v*)^, 
the  movement  is  largely  due  to  the  lower- frequency  components. 
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equations  for  compressible  flow  must  be  reconciled  with  the  incompress¬ 
ible  case  when  the  density  remains  constant.  If  one  applies  a  similarity 
consideration  in  a  very  naive  manner,  the  apparent  conclusion  is  that 


e  %  —  Si  0(v/l) 

Sx 


where  V  and  L  are  the  similarity  scales.  It  would  then  be  difficult 
to  produce  a  consistent  theory  aside  from  making  the  postulate  that 


e  =  0  if  p  =  Constant 
e  =  0(v/l)  if  p  ^  Constant 

To  say  the  least,  such  a  postulate  would  he  too  arbitrary.  A  more 
satisfactory  concept  is  to  regard  e  as  an  entity  by  itself.  Although 
composed  of  terms  like  5u'/Sx,  the  net  effect  of  all  three  terms  com¬ 
bined  together  cause  e  to  behave  in  its  own  way.  The  equation  of 
continuity  (15)  tells  how  this  entity  e  is  related  to  the  density 
variation. 


One  is  now  in  a  position  to  write  down  the  differential  equations 
for  the  lower -frequency  fluctuating  components.  From  equations  (ll^-)  to 
(16),  with  partial  linearization  and  separation  of  lower-  and  higher- 
frequency  components,  the  following  are  obtained: 


^(„  *  .  -i  |_(p .  p,.) 


p  Bx 


v^'  =  -i  — (p  +  p?') 
Dt  ^  pdy^ 


p  Sz 


(21) 


5t(p  ^  Pz')  ^  ° 


(22) 
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^(t  +  T2')  +  (p  +  P2')e 


1J.€  +  (1^6)7' 


Here  —  is  to  be  understood  as 
Dt 


D  5  /-  ,\  S  ,5  .  S 


From  the  equation  of  state  (l7)^  one  gets 


p  =  RpT 


V  ^  V' 

p  p  f 


where  the  double  correlation  p’T*  is  omitted  since 


pT  p  ,  T  ^ 


In  the  next  section^  similarity  scales  will  be  derived  from  the  above 
differential  equation  in  parallel  with  the  incompressible  theory.  The 
mean  velocity  and  temperature  distributions  will  then  follow. 


SIMILARITY  SCALES 


Based  on  the  system  of  differential  equations  (21)  to  (23)  for  the 
lower-frequency  components  of  the  fluctuations,  it  is  now  possible  to 
deduce  the  similarity  scales  for  these  lower- frequency  components  and 
build  upon  them  a  theory  of  the  compressible  boundary  layer  in  parallel 
with  the  incompressible  case.  For  the  sake  of  simplicity,  let  the  sub¬ 
script  I  be  dropped  and  rewrite  the  equations  as:  ' 


D  /- 


(ui  +  Ui') 


P 


^  +  p’),  1  =  1,  2,  3 
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Cp  ^  (T  +  T') 
^  Dt 


R  D_ 
p  Dt 


P’ 


=  V€  +  (ve) ' 


(28) 


where 


p  =  RpT 
£l  =  £l  + 

P  P  T 


(29) 


e 


and 


Dt 


I  \  S  ,  S 
^  Sx  ^  ^ 


+  v’ 


bz 


Equation  (28)  is  obtained  by  eliminating  e  in  the  energy  equation  by 
means  of  the  continuity  relation. 

For  the  case  of  boundary  layer  over  a  flat  plate,  consider  the  state 
of  affairs  at  large  distances  from  the  leading  edge.  There,  if  the  slope 
of  the  outer  edge  of  the  boundary  layer  be  assumed  small  as  is  usually 
done  for  the  incompressible  case, 3  a  first  approximation  will  be  to  regard 


3ln  the  incompressible  case,  a  popular  formula,  though  not  too  well 
substantiated,  for  the  growth  of  the  boundary -layer  thickness  is 

6  «  xV5 

hence  - - >0  as  x- — >00.  See,  e.g.,  reference  21. 
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the  mean  flov  within  the  boundary  layer  as  a  parallel  flow  varying  only 
in  the  normal  direction,  namely 

u  =  u(y) 

V  =  w  =  0 

The  mean  pressure  is,  of  course,  constant  and  eqiml  to  the  value  in  the 
main  stream.  In  considering  local  similarity,  the  observer  may  be  taken 
to  move  with  the  local  mean  motion,  again  as  in  the  incompressible  case. 
Thus, 


u 


1  2  d^ 

2  ^  dy2 


+ 


D_ 

Dt 


du  ,  1  2 

y  —  +  —  y 
dy  2^ 


+  v' 


(30) 


and  so  on;  the  right-hand  side  of  equation  (26)  becomes 

_1  ^ 

P 

and  the  dilatation  reduces  to 

Buj  * 


One  may  recall  that,  at  this  point,  for  the  incompressible  case 
the  procedure  was  to  use  the  vorticity  equation  by  eliminating  the 
pressure  fluctuation  terms  on  the  right-hand  side  of  equation  (26)  from 
each  other.  Such  a  step  is  not  readily  useful  here  because  the  mean 
density  now  varies  with  y.  Instead,  by  cross -multiplication,  there 
follows 


Sp"  Dv'  5p*  D  /-  ,  ^ 

- -  - (u  +  u*) 

ox  Dt  Sy  Dt 


(31) 


and  so  forth.  Hence,  if  a  scale  can  be  found  for  the  pressure  fluctu- 
equation  (31)  indicates  that  the  velocity  and  time  scales  are 
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determined  by  a  relation  identical  with  that  for  the  incompressible  case 
In  other  words,  these  scales  should  come  from  a  consideration  of  the 

operator  D/Dt.  Expanding  ^  (u  +  uV), 


-(u*u')  .  V  - 


'  du  ,  1  2  ^  .  I  1  .  t  .  I  Su* 

7  ^  2  ^  dy2  ""  j  ax  ^  ay  az 


Let  the  fluctuating  motion  have  a  length  scale  Iq,  a  time  scale  Iq, 
and  a  velocity  scale  v^,  where 


^o  “  ^o/'to 


After  substitution, 


L.  /- 

Dt 


d^u 

dy2 


+  . 


vi'  + 


y2  ^o  d^u  ^  ^  \ 

+  .  .  .  +  VqVi'  - 

Iq  dy2  /  2( 


Vo 

2o  ati 

-  -  f-  ^0  -  - 

\2o  °  <iy 

5Un  * 

1 

aui*'' 

^oWl 

azi  ^ 

where  subscript  1  denotes  the  normalized  quantity: 

(UlS  vi')  =  (u',  v’)/vo 
(Xly  yiy  Zl)  =  (x,  y,  z)/Zo 


The  similarity  solutions  require 


^i'  =  ^i'(xi^  yi^  21,  ti) 
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and  so  forth.  Therefore^  the  mean  flow  quantities  must  be  separable 
from  the  normalized  fluctuating  quantities  in  the  above  expression. 
Exactly  as  in  the  incompressible  case,  it  follows  that 


j 


(32) 


and,  if  the  mean  velocity  distribution  is  of  predominant  importance,  so 
that  the  term  involving  must  be  retained. 


du  /d^ 

<iy/ 


(33) 


Applying  the  same  procedure  to  Dv'/Dt,  the  same  set  of  scales  is 
obtained.  Equation  (31)  is  seen  to  be  satisfied  by  such  a  choice  of 
scales,  provided  a  scale  for  pressure  fluctuation  exists. 

The  pressure  scale  is  easily  deduced  from  equation  (26).  Let  the 
scale  be  n.  One  has 


1  ^ 

P  Iq 


—  (u  +  u' ) 
Dt 


Vq  Dui' 


to  Dti 


with  Duj^'/Dt]^  representing  the  normalized  expression  of  ^  (u  +  u'). 
Thus, 


«  «  pvq^  (3^) 

As  the  next  step,  the  energy  equation  (28)  may  be  considered  to 
establish  a  possible  scale  for  the  temperature  fluctuation  T^ .  A  naive 
approach  of  again  assuming  a  single  scale  for  T’  and  substituting  the 
rest  by  the  results  (32)  to  (3^)  will  lead  to  additional  restrictions  on 
the  mean  flow  and  rule  out  a  self-consistent  description.  This  approach 
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must  therefore  be  abandoned.  In  view  of  the  linearity  of  the 
variable  T'  in  equation  (28),  however,  one  can  at  least  formally  split 
the  temperature  fluctuation  T’  into  three  parts  and  rewrite  equa¬ 
tion  (28)  into  a  system  of  equations: 

’0  (35) 


T 


R  D 


P  Dt  2  p  Dt 


p'  =  0 


(36) 


D 

Dt 


V€  +  (vc) ' 


(37) 


Equation  (35)  is  the  same  as  the  usual  incompressible  heat-transfer  equa- 
tion  in  turbulent  flow.  Consequently,  may  be  interpreted  as  the 

temperature  fluctuation  caused  by  the  mixing  of  fluid  elements  from 
strata  of  different  temperatures,  due  to  the  turbulent  velocity  fluctu¬ 
ations.  Alternatively,  one  may  regard  equation  (35)  as  describing  that 
a  certain  portion  of  the  temperature  is  attached  to  each  fluid  element 
during  the  process  of  mixing.  Equation  (36)  describes  the  contribution 
to  temperat'ure  fluctuation  by  the  pressure  field  acting  on  the  volume 
change  of  the  fluid  elements  -  the  compressibility  effect.  Such  an 
interpretation  is  most  clearly  seen  when  the  pressure  term  is  converted 
back  to  the  form  of  equation  (35)  involving  explicitly  the  dilatation  e. 
The  third  part  in  equation  (37)  obviously  represents  the  heating 

effect  of  the  viscous  dissipation. 


One  may  now  introduce  different  scales  63^,  @2^  ®3 

fluctuations  Tp',  T2S  ^3'*  Expansion  of  equation  (35)  leads  to 


,27?r 


,/dT  ^  y  ,  d  T  ^ 

V^Vn  '  I  -  +  - :r  +  .  . 

dy 


■  I,  dy2 


Hence, 


/  *0 

stW;  •  • 

dT 

°  dy 

=  0 


(38) 
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and,  if  the  higher  derivatives  of  the  mean  temperature  distribution  are 
of  predominant  importance. 


In  a  similar  way,  by  expanding  equations  (36)  and  (37)^ 

“  ■'^o^ 


(39) 


(40) 


a  to[v€  +  (ve)i] 


(41) 


It  was  previously  shown  that,  for  the  dissipation,  both  the  mean  and  the 


fluctuating  parts  are  related  to  terms  like 


/ 


and  are  of  the  same 


order  of . magnitude .  Taylor’s  expression  of  balancing  the  energy  yields 
now 


ve  « 


Substitution  into  equation  (^1)  shows  that  0^  is  oi*  same  form  as 

02  •  There  are  but  two  scales  for  the  temperature  fluctuation  given  by 
equations  (38)  and  (4o). 


The  scales  of  the  density  fluctuation  must  satisfy  the  second 
relation  of  equations  (29).  Again  one  has  a  linear  form  and  may  split 
p’  into  several  parts: 


Denoting  the  scales  for  p3_’,  P2',  •  •  •  ^i^  1*2,  .  .  one 

recasts  the  above  into  a  system  of  equations: 


^1 


oc  ^  0- 


oc 


<ay 


(42) 
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^2  '  I 


m 


Thus  T2  and  also  are  proportional.  With  the  help  of  the  first 

relation  of  equations  (29)>  equation  (^2)  may  be  rewritten  as 


r 


1 


a 


(45) 


Equations  (33)  and  (4-5)  are  the  expressions  for  the  two  scales  of  p'. 

It  is  desirable  at  this  moment  to  bring  attention  to  the  entity 
"dilatation."  The  equation  of  continuity  (27)  has  not  been  used  in  the 
previous  derivation  of  the  scales.  One  could  formally  write  e  oc  Vq/Zq 
and  find  that  the  above  scales  satisfy  equation  (27)  as  well.  A  proce¬ 
dure  like  this,  however,  is  not  justifiable  on  the  ground  of  the  argu¬ 
ments  in  the  section  "Extension  of  Similarity  Theory  to  Compressible 
Turbulent  Boundary  Layer." 

All  the  scales  of  the  lower -frequency  fluctuations  have  now  been 
related  to  the  mean  flow.  There  remains  one  point  which  needs  some 
clarification.  The  question  is:  For  the  length  scale  Zq,  which  of  the 

expressions,  (33)  or  (39),  should  be  adopted?  Evidently,  if  they  are 
compatible,  there  must  be 
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Integrating, 


On  the  other  hand,  it  will  be  seen  in  reference  7  that  a  relation 

between  u  and  T  may  be  obtained  from  averaging  the  equations  of 
motion  and  energy .  The  relation  is  further  independent  of  the  choice 
of  Zq.  If  the  expressions  of  Zq  from  equations  (33)  and  (39)  are 

equivalent,  the  relation  must  also  be  reconcilable  with  equation  (^7). 

Such  is  the  case  for  incompressible  flow  when  Reynolds'  analogy  is  valid 
(i.e.,  ai/a2  =  1).  But,  in  general,  one  can  hardly  expect  the  com¬ 
patibility  relation  (^7)  to  be  fulfilled,  except  perhaps  approximately. 

It  will  be  seen  later  that  the  choice  of  the  form  of  Zq  determines  the 

distribution  in  the  space  coordinate  (y-direction) .  At  present,  one  will 
be  contented  with  the  notion  that  the  choice  of  equation  (33)^  say,  means 
more  emphasis  on  the  satisfaction  of  the  momentum  relations  at  the  expense 
of  the  energy  relations,  and  vice  versa  for  the  choice  of  equation  (39). 
Some  additional  discussions  are  included  in  the  following  section. 


DISCUSSION 


Before  a  detailed  treatment  of  the  theory  of  turbulent  boundary 
layer,  it  may  be  worth  while  to  remind  oneself  of  all  the  approximations 
involved  and,  hence,  the  limitations  to  be  Introduced  in  the  investi¬ 
gation. 

The  preceding  theory  for  the  turbulent  boundary  layer  involves  mainly 
the  following  approximations: 

(1)  The  turbulence  pattern:  The  turbulence  pattern  is  idealized  to 
be  such  that  the  small  (high-frequency)  dissipative  eddies  are  separable 
from  the  larger  (lower -frequency)  energy-carrying  eddies.  In  accordance 
with  Kolmogoroff  (reference  19),  the  small  eddies  are  isotropic  and 
participate  in  an  equilibrium  governed  by  the  kinematic  viscosity  v 
and  the  rate  of  dissipation  e.  The  larger  eddies  are  assumed  also  to 
participate  in  an  equilibrium,  governed  by  the  local  conditions.  Here 
one  sees  that  three  approximations  are  involved: 

(a)  The  separability  of  small  and  larger  eddies 
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(Id)  The  possibility  of  similarity  for  the  larger  eddies 

(c)  The  "local^’  character  of  the  similarity  for  the  larger  eddies 

In  order  that  the  small  dissipative  eddies  may  be  separated^  an  essential 
condition  is  that  the  turbulence  Reynolds  number  be  high.  This^  in 
general,  will  follow  if  the  mean-flow  Reynolds  number  is  high,  for  it  is 
a  commonly  accepted  notion  that  a  higher  mean-flow  Reynolds  number  creates 
more  small  eddies.  In  this  respect,  a  higher  main-stream  velocity  out¬ 
side  the  boundary  layer  should  check  better  with  the  theory.  The  simi¬ 
larity  for  the  larger  eddies  is  mainly  an  intuitive  hypothesis  and  cer¬ 
tainly  needs  eD^erimental  verification.  The  arguments  for  its  plausi¬ 
bility  are  given  in  reference  6.  In  the  very  simple  configuration  of  a 
flat-plate  boundary  layer  without  pressure  gradient,  the  chances  for  the 
realization  of  similarity  are  certainly  on  the  favorable  side.  One  may 
even  take  a  heuristic  point  of  view  and  regard  the  whole  theory  as  being 
an  attempt  to  seek  a  similarity  solution,  if  at  all  possible.  The 
behavior  of  the  small  eddies  is  decidedly  different,  but  the  order-of- 
magnitude  study  in  reference  6  indicates  that  the  original  formulation 
by  Von  Karman  in  1930  (reference  22)  remains  valid  if  the  fluctuations 
are  interpreted  as  the  lower -frequency  ones  only.  The  scales  for  length 
and  time  are  determined  from  the  local  state  of  the  mean  flow  by 
restricting  attention  to  the  immediate  neighborhood  of  the  point  in 
question.  It  is  well-known  that  the  length  scale  from  such  considerations 
is  only  moderately  small  in  comparison  with  the  boundary -layer  thickness. 
Experiments  in  the  National  Bureau  of  Standards  for  boundary  layer  with 
a  pressure  gradient  fxirther  showed  the  existence  of  significant  large 
eddies  (reference  23).  However,  there  is  the  consolation  that  for  a 
class  of  the  distribution  profile,  pointed  out  by  Prandtl  (reference  21, 
p.  132)> 


u  =  A(y  +  +  C  (W) 

\ 

where  A,  B,  and  C  are  constants,  the  magnitude  of  the  length  scale 
is  of  no  consequence.  This  is  true  because  the  higher  derivatives,  which 
bring  in  the  conditions  away  from  the  point,  will  form  length  scales 

proportional  to  the  one  determined  by  du/dy  and  d^/dy^.  The  loga¬ 
rithmic  distribution  belongs  to  the  class  (i^-8)  and  so  do  the  usual  power 
(1/7,  1/9^  etc.)  laws. 

(2)  The  boundary -layer  growth:  The  growth  of  the  thickness  of  the 
boundary  layer  has  been  neglected  in  the  present  theory,  as  everything 
is  assumed  to  depend  on  the  coordinate  y  only.  Such  a  neglect  is 
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equivalent  to  the  stipulation  that  dS/dx  «  1,  &  being  the  boundary- 

layer  thickness.  In  the  incompressible  case,  semi  empirically  the  gen¬ 
erally  accepted  formula  is 

5  oc  xV5 

which  shows  that,  although  the  rate  of  decrease  of  the  slope  is  rather 
slow,  db/dx  could  be  as  small  as  one  wishes  at  distances  very  far  from 
the  leading  edge.  In  applying  the  theory  to  moderate  distances  from  the 
leadiing  edge,  one  must  be  aware  of  this  limitation  when  compared  with 
experiments . 

(3)  The  two  choices  of  the  length  scale:  With  the  extension  to 
the  compressible  case,  two  possible  choices  of  the  length  scale  are 
found,  depending  on  whether  the  equations  of  momentum  or  the  equation 
of  energy  is  to  be  better  satisfied.  Unless  the  two  scales  are  essen¬ 
tially  identical,  in  the  manner  discussed  at  the  end  of  the  section 
"Similarity  Scales,"  an  added  degree  of  approximation  is  involved, 
because  one  of  the  expansions  has  to  be  cut  down  to  the  order  of  the 
first  derivative  only.  In  this  situation,  the  validity  therefore  hinges 
on  whether  the  second  derivative  i£  sufficiently  small  in  comparison 
with  the  first  for  either  u  or  T.  In  other  words,  the  approximation 
will  still  be  good  if  one  of  the  two  scales  is  much  larger  than  the 
other,  for  example. 


du  /d^ 

cly/dy2 


«  — 


dT  /d^T 


dy/  dy^ 


so  that,  for  local  behavior,  only  the  smaller  scale  (here 

\  ^y/<3y2/ 

controls  the  motion.  As  an  extreme  case  for  the  purpose  of  illustration, 
suppose  that  a  large  amount  of  heat  were  taken  away  at  the  flat  plate  to 
reduce  the  wall  temperature  to  that  in  the  free  stream.  The  temperature 
profile_in  the  boundary  layer  would  be  approximately  uniform.  The 
term  dT/dy,  and  so  forth  enters  the  determination  of  the  length  scale 
through  the  convective  action  and  is  obviously  unimportant  in  this  case. 
One  must  then  choose  the  length  scale  from  the  velocity  profile.  In  the 
general  case,  it  can  only  be  hoped  that  the  actual  distributions  do  not 
depart  too  much  from  the  relation  (46),  which  holds  for  approximately 
identical  scales. 


(4)  The  effect  of  pressure  gradient:  It  is  interesting  to  observe 
that  the  Introduction  of  a  pressure  gradient  renders  impossible  the 
extension  to  the  compressible  case.  This  is  so  because  the  scale  for 
the  pressure  fluctuation  can  no  longer  be  determined  as  related  only  to 
the  velocity  scale  (cf.  equation  (34)).  Without  a  similarity  among  the 
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pressure  fluctuations^  there  follows  no  similarity  for  the  temperature 
and  the  density  fluctuations.  The  entire  scheme  breaks  down.  Difficulty 
like  this  does  not  influence  the  incompressible  case,  where  the  pressure 
fluctuation  is  of  no  consequence  and  of  no  interest.  A  remark  may  be  in 
order  about  Von  Karman’s  initial  derivation  in  considering  a  channel 
flow.  His  derivation  started  from  the  vorticity  equation.  Such  a  pro¬ 
cedure  is  still  permissible  in  the  present  theory  when  p  remains  con¬ 
stant,  so  that  pressure  terms  can  be  eliminated  by  cross -differentiation 
of  equation  (26).  The  same  similarity  theory  would  then  be  obtained. 

When  the  results  of  the  present  theory  for  a  flat-plate  boundary 
layer  are  to  be  applied  to  the  case  of  an  airfoil,  the  additional  restric 
tion  to  the  possible  degree  of  agreement  should  be  kept  in  mind. 


Massachusetts  Institute  of  Technology 

Cambridge,  Mass.,  December  27,  1950 
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SYMBOLS 


^1^  ^2 


e 


k 

I 


P 

P^Pl 


constants 

specific  heat  at  constant  pressure  and  constant  volume, 
respectively 

coefficient  of  heat  conductivity  of  fluid 
mixture  len^h  or  similarity  scale  of  length 
similarity  scale  of  length 

pressure 

mean  pressure  and  mean  pressure  in  free  stream,  respectively 


ratios  of  various  heat  conduction  and  viscous  dissipation 
terms  to  tiirbulent  transfer  term  in  equation  of  energy, 
discussed  in  section  "Energy  Equation  and  Dissipation 
Term"  *  . 

similarity  scales  of  three  components  of  density  fluctu¬ 
ation,  defined  by  equations  (42)  to  (44) 

t  time 

similarity  scale  of  time  for  fluctuations 

u,v,w  velocity  components  in  x-,  y-,  and  z-directions, 

respectively 

x,y,z  Cartesian  coordinates;  x,  axis  in  direction  of  plate  and 

free  stream;  y,  axis  normal  to  plate;  and  z,  axis 
parallel  to  leading  edge  of  plate 

L 


characteristic  lexigth 
characteristic  Mach  niunber 
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3h 


R- 


R 


uv 


R. 


T, 


m 


AT. 


m 


U 


U 


m 


€,€ 


X 

V 

Tt 

P 

T 

d) 


gas  constant  in  equation  of  state 
characteristic  Reynolds  number 

correlation  constant  between  u’  and  v’ 

turbulence  Reynolds  number 

temperature 

characteristic  mean  temperature 

characteristic  mean  temperature  difference 

mean  temperature  at  wall 

mean  velocity  component  in  x-direction 
characteristic  mean  velocity 

ratio  of  specific  heats 
thickness  of  boundary  layer 

rate  of  dissipation  and  its  fluctuating  component^ 
respectively 

similarity  scales  of  temperature  fluctuation  defined  by 
equations  (38),  (to),  and  (4l),  respectively 

Taylor's  microscale  of  turbulence 

coefficient  of  viscosity 

coefficient  of  kinematic  viscosity 

similarity  scale  of  pressure  fluctuation 

density  of  fluid 

shearing  stress 

dissipation  function  in  equation  of  energy 
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Subscripts: 

I  low-frequency  part  of  fluctuations 

h  high-frequency  part  of  fluctuations 

1  normalized  fluctuating  components^  for  example^ 

f^l'  =  M-'/vo7  x'  =  x/Zq,  and  so  forth;  also,  quantities 
in  free  stream 

Barred  quantities  always  represent  mean  values;  primed  quantities 
represent  fluctuations. 
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Von  Karman’s  theory  in  incompressible  flow,  simi-  Von  Karman’s  theory  in  incompressible  flow,  simi¬ 
larity  scales  for  all  the  flow  variables  are  derived.  larity  scales  for  all  the  flow  variables  are  derived. 

Two  possible  length  scales  are  foimd  and  the  signifi-  Two  possible  length  scales  are  found  and  the  signifi¬ 
cance  of  this  possibility  is  discussed.  cance  of  this  possibility  is  discussed. 


•d  3 

If 

U  V} 
d) 

'I 

-I 

w  ^ 


W  O  o 

O*  4) 


